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The cross sections of е e^ — тл he at center-of-mass energies from 3.896 to 4.600 GeV are mea- 
sured using data samples collected with the BESIII detector operating at the Beijing Electron Positron Col- 
lider. The cross sections are found to be of the same order of magnitude as those of e^ e^ — ata J/w and 
ete — тїт (2S), but the line shape is inconsistent with the Y states observed in the latter two modes. 
Two structures are observed in ће e* e^ — тл he cross sections around 4.22 and 4.39 GeV/c’, which we 
call Y (4220) and Y (4390), respectively. A fit with a coherent sum of two Breit-Wigner functions results in a 
mass of (4218.4* 5? + 0.9) MeV/c? and a width of (66.0* 12? + 0.4) MeV for the Y (4220), and a mass of 
(4391.5* 83 + 1.0) MeV/c? and a width of (139.5*362 + 0.6) MeV for the Y (4390), where the first uncer- 
tainties are statistical and the second ones systematic. The statistical significance of Y (4220) and Y (4390) is 


10c over one structure assumption. 


PACS numbers: 14.40.Rt, 14.40.Pq, 13.66.Bc, 13.25.Gv 


In the last decade, a series of charmonium-like states 
have been observed at e*e~ colliders. These states chal- 
lenge the understanding of charmonium spectroscopy as well 
as QCD calculations а. According to potential mod- 
els, there are five vector charmonium states between the 1D 
state (3770) and 4.7 GeV/c?, namely the 35, 2D, 45, 
3D, and 55 states Ш. From experimental studies, besides 
the three well-established structures observed in the inclu- 
sive hadronic cross section [B], i.e., (4040), (4160), and 
(4415), five Y states, i.e., Y (4008), Y (4230), Y (4260), 
Y (4360), and Y (4660) have been reported in initial state ra- 
diation (ISR) processes e*e^ > «rnm ^x J/i оге е — 
звт їл (2,6) at the B factories or in the di- 
rect production processes at the CLEO and BESIII experi- 
ments The overpopulation of structures in this re- 
gion and the mismatch of the properties between the poten- 
tial model prediction and experimental measurements make 
them good candidates for exotic states. Various scenarios have 
been proposed, which interpret one or some of them as hybrid 
states, tetraquark states, or molecular states [14]. 


The study of charmoniumlike states in different produc- 
tion processes supplies useful information on their proper- 
ties. The process єтє — тл Һе was first studied by 
CLEO [15] at center-of-mass (c.m.) energies \/s from 4.000 
to 4.260 GeV. A 10c signal at 4.170 GeV and a hint of a rising 
cross section at 4.260 GeV were observed. Using data sam- 


ples taken at 13 c.m. energies from 3.900 to 4.420 GeV [1б], 
BESIII reported the measurement of the cross section of 
e*e- — «*z-h, [17]. Unlike the line shape of the pro- 
cess ete — mm J/y, there is a broad structure in the 
high energy region with a possible local maximum at around 
4.23 GeV in e*e- — m«*m«-h,. Based on the CLEO mea- 
surement at \/s = 4.170 GeV and the BESIII measurement, 
two assumptions were made to describe the cross section in 
Ref. [18]. In both assumptions, a narrow structure exists at 
around 4.23 GeV, while the situation in the high energy re- 
gion is unclear due to the lack of experimental data. 


In this Letter, we present a follow-up study of етет — 
тїт he at c.m. energies from 3.896 to 4.600 GeV us- 
ing data samples taken at 79 energy points with the 
BESIII detector (20). Two resonant structures are observed 
at 4/8 = 4.22 and 4.39 GeV [hereafter referred to as Y (4220) 
and Y(4390)]. The integrated luminosity at each energy 
point is measured with an uncertainty of 1.0% using large- 
angle Bhabha events рд. There аге 17 energy points 
where the integrated luminosities are larger than 40 pb! (re- 
ferred to as "XY Z data sample" hereafter), while the inte- 
grated luminosities for the other energy points are smaller 
than 20 pb~! (referred to as “R-scan data sample" hereafter). 
The c.m. energies for the XY Z data sample are measured 
with етет + ^trsn/rsnH* HT events with an uncertainty of 
+0.8 MeV [23], which is dominated by the systematic uncer- 
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tainty. A similar method is used for the R-scan data sample 
with multihadron final states (24h. 


In this study, the he is reconstructed via its electric- 
dipole transition he — yne with ne —› Xi, where 
X; is one of 16 exclusive hadronic final states: рр, 
(пт), 2(KtK-), ntn KK, mm рр, 3x"), 
2m m KK, KOK at, KOK^mm a, KOK a, 
ppr®, K*K m, ntan, 2(n*'n )n, път nr, and 
2(r*m- m?). Here, the А0 is reconstructed using its decay 


to rtr, and the 7° and 7 from the yy final state. 


Monte Carlo (MC) simulated events are used to optimize 
the selection criteria, determine the detection efficiency, and 
estimate the possible backgrounds. The simulation of the 
BESIII detector is based on GEANT4 and includes the 
geometric description of the BESIII detector and the detec- 
tor response. For the signal process, we use an MC sample 
for ete — петт h, process generated according to phase 
space. ISR is simulated with KKMC with a maximum en- 
ergy for the ISR photon corresponding to the тл he mass 
threshold. 


We select signal candidates with the same method as that 
described in Ref. [17]. Figure[I] shows the scatter plot of the 
invariant mass of the 7. candidate versus the one of the he 
candidate and the invariant mass distribution of уле in the 
Ne signal region for the data sample at ys = 4.416 GeV. 
A clear he — Yne signal is observed. The ne signal re- 
gion is defined by a mass window around the nominal ne 
mass B), which is within +50 MeV/c? with efficiency about 
84% (+45 MeV/c? with efficiency about 80%) from MC sim- 
ulation for final states with only charged ог КО particles (for 
those including 7° or 7). 


We determine the number of път he signal events (поре) 
from the yne invariant mass distribution. For the XY Z data 
sample, the yne mass spectrum is fitted with the MC simu- 
lated signal shape convolved with a Gaussian function to re- 
flect the mass resolution difference between the data and MC 
simulation, together with a linear background. The fit to the 
data sample at 4/5 = 4.416 GeV is shown in Fig.[I] The tail 
on the high mass side is due to events with ISR (ISR photon 
undetected); this is simulated with KKMC in MC simulation, 
and its fraction is fixed in the fit. For the data samples with 
large statistics (ys = 4.226, 4.258, 4.358, and 4.416 GeV), 
the fit is applied to the 16 т decay modes simultaneously with 
the number of signal events in each decay mode constrained 
by the corresponding branching fraction (27). For the data 
samples at the other energy points, we fit the mass spectrum 
summed over all п. decay modes. For the R-scan data sample, 
the number of signal events is calculated by counting the en- 
tries in the ^, signal region [3.515, 3.535] GeV/c? (155) and 
the entries in the A, sideband regions [3.475, 3.495] GeV/ c? 
and [3.555, 3.575] GeV/c? (п) using the formula nu = 
nis — fnde, Here, the scale factor f = 0.5 is the ratio of the 
size of the signal region and the background region, and the 
background is assumed to be distributed linearly in the region 


of interest. 
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FIG. 1. The Myne distribution in the 7. signal region of 4.416 GeV 
data. Points with error bars are the data and the curves are the best fit 
described in the text. The inset is the scatter plot of the mass of the 
Ne candidate Mn, versus the mass of the he candidate Myn, for the 
same data sample. 


The Born cross section is calculated from 


obs 
B. np. 
L(1 + 6)|1 + 128, У) e; Boi) 
where mere is the number of observed signal events, £ is 


e 


the integrated luminosity, (1 + 8) is the ISR correction factor, 
|1 + П? is the correction factor for vacuum polarization (23h, 
Ві is the branching fraction of he — yne Ві, є; and B»(1) are 
the detection efficiency and branching fraction for the 7 th ne 
decay mode (27h, respectively. The ISR correction factor is 
obtained using the QED calculation as described in Ref. 
and taking the formula used to fit the cross section measured 
in this analysis after two iterations as input. The Born cross 
sections are summarized in the Supplemental Material 
together with all numbers used in the calculation of the Born 
cross sections. The dressed cross sections (including vacuum 
polarization effects) are shown in Fig.[2]with dots and squares 
for the R-scan and X Y Z data sample, respectively. The cross 
sections are of the same order of magnitude as those of the 
ete” > ata J/wandete” > ntn (2S) [412], but 
follow a different line shape. The cross section drops in the 
high energy region, but more slowly than for the e*e^ — 
vm J/w process. 


Systematic uncertainties in the cross section measurement 
mainly come from the luminosity measurement, the branching 
fraction of he — уђе, and ne — Х;, the detection efficiency, 
the ISR correction factor, and the fit. The uncertainty due to 
the vacuum polarization is negligible. The uncertainty in the 
integrated luminosity is 1% at each energy point. The uncer- 
tainty sources for the detection efficiency include systematic 
uncertainties in tracking efficiency (196 per track), photon re- 
construction (1% per photon), and I$ reconstruction (1.2% 
per К 2); Further uncertainties arise from the 70/7) mass win- 
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dow requirement (1% per 70/2), the uis requirement, ne pa- 
rameters, and line shape, possible intermediate states in the 
7h, and ттлт mass spectra, intermediate states in ne de- 
cays (included in the uncertainty from the branching fraction 
of ne — X;), and the limited statistics of the MC simulation. 


The uncertainty due to the Mid requirement is estimated 
by correcting the helix parameters of the simulated charged 
tracks to match the resolution found in data, and repeating the 
analysis Bd). Uncertainties due to the т. parameters and line 
shape are estimated by varying them in the MC simulation. 
When producing MC events for the e*e^ — т*л Ре pro- 
cess through the intermediate states Z.(3900) or Z.(4020), 
the parameters of the 2.(3900) and 2„(4020) are fixed to the 
average values from the published measurements 


[33]. The quantum numbers of both Z.(3900) and Z.(4020) 


are assumed to be J” = 1+. The differences in the efficiency 
obtained from phase space MC samples and those with inter- 
mediate Ze states are taken as the uncertainties from possi- 
ble intermediate states in the 7*h, system. The uncertainty 
from intermediate states in the тт system is estimated by 
reweighting the тл mass distribution in the phase space 
MC sample according to the data, and the resulting difference 
in the efficiency is considered as uncertainty. The uncertain- 
ties due to data and MC differences in the detection efficiency 
are determined to be between 5.596 and 10.896, depending on 
the т. decay modes and the c.m. energy. Combining the un- 
certainties for the branching fractions of 7. decays bà. the 
uncertainties for the average efficiency SE є;182(1) are be- 
tween 6.4% and 9.1% depending on the c.m. energy. 


The uncertainty in the ISR correction is estimated as de- 
scribed in Ref. 1l. Uncertainties due to the choice of the 
signal shape, the background shape, the mass resolution, and 
fit range are estimated by changing the he and ne resonant 
parameters and line shapes in the MC simulation, changing 
the background function from a linear to a second-order poly- 
nomial, changing the mass resolution difference between the 
data and the MC simulation by 1 standard deviation, and by 
extending or shrinking the fit range. 


Assuming all of the sources are independent, the total sys- 
tematic uncertainty in the л+л he cross section measurement 
is determined to be 9.4%-13.6% depending on the c.m. en- 
ergy. The uncertainty in B; is 11.8% [all common to all en- 
ergy points, and quoted separately in the cross section mea- 
surement. Altogether, the quadratic sum of the common sys- 
tematic errors at each energy point accounts for about 95% of 
the total systematic error. 


A maximum likelihood method is used to fit the dressed 
cross sections to determine the parameters of the resonant 
structures. The likelihood is constructed taking the fluctua- 
tions of the number of signal and background events into ac- 
count (the definition is described in the Supplemental Mate- 
rial ). Assuming that the тту signal comes from two 
resonances, the cross section is parametrized as the coherent 
sum of two constant width relativistic Breit- Wigner functions, 


e(m) =|Bi(m)- 


P(m) i P(m) 2 
PL) + e? Bo(m) - V Pun! 


where B;(m) = 127(г..В);Г, /(m? = м? + іМ;Г;) 
with j = 1 or 2 is ће Breit-Wigner function, and P(m) 
is the three-body phase space factor. The masses M;, the 
total widths Гу, the products of the electronic partial width 
and the branching fraction to път he (Г eeD);, and the rel- 
ative phase ф between the two Breit-Wigner functions are 
free parameters in the fit. Only the statistical uncertainty 
is considered in the fit. There are two solutions from the 
fit, one of them is shown in Fig. The second solu- 
tion is very close to the one shown here. This can been 
proved analytically using Eq.(9) in Ref. (34), which relates 
the two solutions from the fit when a sum of two coher- 
ent Breit-Wigner functions is used. The parameters deter- 
mined from the fit are M; = (4218.47?) MeV/c?, Ty = 
(66.01 53°) MeV, and (Г.г): = (4.6129) eV for Y (4220), 
Мз = (4391.5*63) MeV/c?, T» = (139.5*192) MeV, and 
(Г..8)ә = (11.6159) eV for Y (4390). The relative phase 
$ is (З) rad. The correlation matrix of the fit parame- 
ters shows large correlation between the (Г..8); and Ф (see 
Supplemental Material (19). 
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FIG. 2. Fit to the dressed cross section of e^ e^ — ata7he with 
the coherent sum of two Breit-Wigner functions (solid curve). The 
dash (dash-dot) curve shows the contribution from the two structures 
Y (4220) [Y (4390)]. The dots with error bars are the cross sections 
for the R-scan data sample, the squares with error bars are the cross 
sections for the X Y Z data sample. Here the error bars are statistical 
uncertainty only. 


The likelihood contours in the mass and width planes for 
Y (4220) and Y (4390) are shown in Fig.[3] together with the 
positions of Y (4230), Y (4260), Y (4360), and (4415) with 
the parameters taken from the latest PDG average [3]. The 
low-lying resonance from the study of e*e^ — петт 7/4 
at BESIII [35], marked as Y (4260)PES!! in the plot, is also 
compared. Y (4260), Y (4360), and (4415) are located out- 
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side the 3c contours, while Y (4230) and Y (4260)PE5!I are 
overlapped with the 3c contour of Y (4220). 


025 T 025 T T T 
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FIG. 3. The likelihood contours in the mass and width planes for 
Y (4220) (left panel) and Y (4390) (right panel). The filled areas 
are up to 3c likelihood contours and the dots with error bars are the 
locations of Y or w states. The parameters of Y (4260)? P6 are taken 
from the PDG average [3] and Y (4260)®S" from the measurement 
of ete~ — път J/y at ВЕЅШ [35]. 


Fitting the dressed cross section with only one resonance 
yields a worse result, the change of the likelihood value from 
two resonances to one resonance is [A(—2InL) = 113.5]. 
Taking the change in the number of degrees of freedom (4) 
into account, the significance for the assumption of two reso- 
nant structures over the assumption of one resonant structure 
is 10c. The fit with the coherent sum of one Breit-Wigner 
function and a phase space term gives a worse result as well, 
the change of the likelihood value is [A(—21nL) = 66.8]. We 
also fit the cross section with the coherent sum of three Breit- 
Wigner functions, or the coherent sum of two Breit-Wigner 
functions and a phase space term. Both assumptions improve 
the fit quality, but the significances of the third resonance and 
the phase space term are only 2.60 and 2.96, respectively. 


The systematic uncertainties in the resonance parameters 
mainly come from the absolute c.m. energy measurement, 
the c.m. energy spread, and the systematic uncertainty on the 
cross section measurement. The uncertainty from the c.m. en- 
ergy measurement includes the uncertainty of the c.m. energy 
and the assumption made in the measurement for the R-scan 
data sample. Because of the low statistics at each energy point 
in the R-scan data sample, we approximate the difference be- 
tween the requested and the actual c.m. energy by a common 
constant. To assess the systematic uncertainty connected with 
this assumption, we replace the constant by a c.m. energy- 
dependent second-order polynomial. The systematic uncer- 
tainty of the c.m. energy is common for all the energy points 
in the two data samples and will propagate to the mass mea- 
surement (0.8 MeV). The changes on the parameters are taken 
as uncertainty. The uncertainty from c.m. energy spread is es- 
timated by convoluting the fit formula with a Gaussian func- 
tion with a width of 1.6 MeV, which is beam spread, measured 
by the Beam Energy Measurement System . The uncer- 
tainty from the cross section measurement is divided into two 


parts. The first one is uncorrelated among the different c.m. 
energy points and comes mainly from the fit to the yne invari- 
ant mass spectrum to determine the signal yields. The corre- 
sponding uncertainty is estimated by including the uncertainty 
in the fit to the cross section, and taking the differences on 
the parameters as uncertainties. The second part includes all 
the other sources, is common for all data points (14.896), and 
only affects the Г../3 measurement. Table [I] summarizes the 
systematic uncertainty in the resonance parameters. 


In summary, we measure the e*e- — ntn” h, Born cross 
section using data at 79 c.m. energy points from 3.896 to 
4.600 GeV. Assuming ће 7*z~h, events come from two 
resonances, we obtain M — (4218.4153 + 0.9) MeV/c?, 
Г = (66.0722 +0.4) MeV, and (TeeB) = (4.64? 940.8) eV 
for Y(4220), and М = (4391.516 + 1.0) MeV/c?, Г = 
(139.5736°3 + 0.6) MeV, and (T44B) = (11.6150 + 1.9) eV 
for Y (4390), with a relative phase of ф = (3.1705 +0.2) rad. 
The first errors are statistical and the second are systematic. 
The parameters of these structures are different from those of 
Y (4260), Y (4360), and (4415) [3]. The resonance param- 
eters of Y (4220) are consistent with those of the resonance 
observed in e* e^ > wx«o [13]. 


The two resonances observed in e* e^ — ntr" he process 
are located in the mass region between 4.2 and 4.4 GeV/c?, 
where the vector charmonium hybrid states are predicted from 
various QCD calculations 169. The mass of Y (4220) 
is lower than that of Y (4260) observed in the e*e^ — 
пт J/w process. The smaller mass is consistent with some 
of the theoretical calculations for the mass of Y (4260) when 
explaining it as a Dı D molecule (40,41. 
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TABLE I. The systematic uncertainty in the measurement of the resonance parameters. c.m. energy1,2 represent the uncertainty from the 
systematic uncertainty of c.m. energy measurement and the assumption made in the c.m. energy measurement for R-scan data sample, 
respectively. Cross section, оь) represents the uncertainty from the systematic uncertainties of the cross section measurement which are un- 


correlated (common) in each energy point. 


Y Y (4220) Y (4390) did 
М (MeV/c?) T (MeV) (Г..8) (еу)| M(MeV/c) IT (MeV) (Te) (eV) 
c.m. energyi(2) 0.8(0.1) —(0.1) —(0.2) 0.8(0.1) —(0.2) —(0.3  |-(0.D 
c.m. energy spread 0.1 0.3 0.3 0.1 0.1 0.7 0.1 
Cross section, (y) 0.1(—) —(-) 0.2(0.7) 0.6(—) 0.5(—) 0.4(1.7  |0.1(—) 
Total 0.2 
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Appendices 


Cross section of ете — ntr” he 


The number of signal events ns. the luminosity £, the product of the initial state radiation correction factor and average 


i 2 and the Born cross section oP for 
XYZ data sample and R-scan data sample are summarized in Table[IIJand Table [III] The average efficiency is smaller for the 
R-scan data sample than for the X Y Z data sample due to the different methods used in determining the number of signal events. 


efficiency (1 + 5) 35,9 1 €B (Ne — Xi), the vacuum polarization correction factor |1 + П 


TABLE II. e*e^ — mtr” he cross sections from XY Z data sample. The first errors are statistical, and the second ones systematic uncertainty 
except the uncertainty in B(h« — yne), and the third errors are from the uncertainty in B(he — yne). 


Vs (GeV) £(pp 1) т» (1+5), «В(пс  Xi9) +1 a” (pb) 
3.8962 526 —15*1? 0.78 1049  —6.87*51 406408 
4.0076 4820 7.61772 1.42 1.044 2.0779 +0.3+0.2 
4.0855 526 1.61% 1.85 1.052 —3.1185 0.4 0.4 
4.1886 AA атт 1.96 1.056 16.61125 +1.8+2.0 
4.2077 546 23.6+12 1.99 1.057 40.7753 +3.9 4.8 
4.2171 54.1 2531$ 2.03 1.057 4327190 +4.2 5.1 
4.2263 10917 669 + 32 2.07 1.056 552+2.6+6.1- 6.5 
4.2417 556. ЭБЕ. 2.17 1056 39.925 3.84.7 
4.70580 825.7 443 +27 228 1.054 463:2.8 5.12 5.5 
4.3079 449. 43707 2.22 1.052 65.0 158 £6.12 7.7 
4.85583 539.8 357+ 24 2.26 1.051 55.8 + 3.7 + 6.2 + 6.6 
4.3874 55.2 3099.7 2.30 1.051 45.512? E 44 5.4 
4.4156 1073.6 685 + 35 2.28 1.053 521+ 2.7 5.9 6.1 
4.4671 109.9 27.1127 2.38 1055 19.3789 +24 -+2.3 
4.5271 1100 —-13*75 2.39 1.055  —0.9*22 +0.1 0.1 
4.5745 47.7 umo 2.37 1055 15.1112 +21+18 
4.5995 566.9 52.0719 2.38 1.055 7.2 t23 0.7 +0.8 
TABLE Ш: e* e^ — ттт he cross sections in R-scan data sample. The errors 
are statistical only. The systematic uncertainty is 18.0%, and common for all 
energy points. 
TTA) 7 (nh!) 458g „side obs [уту үм о p/p, ү үү 11 LN? Brem | 
v/s (GeV) L (pb !) п®® n тр (14+ 6) 052, «Ble + Х:)(%) |L+ HIP o? (pb) 
40974 7254 0 3 -1513 1.37 1.052 —28.07325 
4.1074 7146 1 0 10723 1.38 1052 “ise. 
4.1174 7648 1 1 0555 1.39 1052 88+1°1 
41274 7207 1 4 1.0473 1.40 1.053 -d85'515 
41374 7.268 4 5 1.573% 1.41 1.0827 Этого 
41224 7.774 1 5 -1.512% 1.42 1052 2599541 
41474 7.662 2 0 2071 1.43 1.053  34.0:228 
44574 7.954 1 1 0.515 1.45 1.054 Өлү 
4.1674 18008 2 5 0.511 1.46 1.054 —3.5*222 
4.1774 7309 1 4 -10:25 1.45 1:055 -175° 8% 
41874 7560 0 3 —1.5t)2 1.45 1056 =—25.553)5 
4194 750 4 4 20? 1.45 1057 “3420505 
41974 7.582 8 3 6.5122 1.45 1057 1094178 
42004 6815 1 4 —1.0+?8 1.46 1057 18075. 
4.2034 74638 2 3 0.5720 1.46 1057 Bates 


Continued on next page 
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TABLE III - continued from previous page 


V/s (GeV) L (pb !) n?8 п пә?  (14-6)Y 5, , Bln — Xi) (%) | -- 12 о? (pb) 
42074 7678 5 5 25138 1.47 1.057 4111871 
42124 7368 0 4 ар 1.48 1056 —324*110 
4.2174 7935 5 2 40125 1.49 1057 62.9+579 
42224 8212. 2 3 05132 1.50 1.057 n5 
4.2274 8193 09 2 8012 1.52 1.056 119.6+$2-4 
4234 8273 8 6 5032 1.53 1056 — 734* 55 
42374 7830 7 6 401 1.54 1.056 61.77 33.3 
4240 8571 7 5 T 1.54 o 63.2755 
42424 8487 7 1 6.5129 1.54 1.05 ба 
4.2454 8.554 10 4 80125 1.55 1055 112.54683 
4.2474 85906 9 6 60725 1.55 1.055 8397925 
4254 8.657 12 5 9.518 1.55 1.054 132.0753 
42574 8880 7 8 302 1.55 1054 40.6152 
42604 8629 6 4 401% 1.55 1053 55.759455 
42674 8548 3 4 101 1.55 1.053 та 
4224 8567 3 5 0.5190 1.56 1.053 т.0+570 
4274 8723 7 4 54. 1.56 1:053 Жа 
42824 8.596 11 7 175155 1.57 1.053 103.6+%82 
42874 9.010 6 11 0.5722 1.57 1.053 667550 
42974 8453 8 7 4511 1.57 1052 63.0195 
43074 859 8 9 3.5111 1.57 1.052 48.2757 
4.3174 9342 8 6 5.011 1.57 1052 63.54 353 
4324 8657 7 з 5529 1.57 0051] 75551 
4334 8700 9 3 1,511 1.58 1051 302052 
43474 8542 7 4 5.0111 1.59 1051 68.81583 
43574 803 8 5 5.511 1.60 1051 79.5 5573 
4.3674 8498 8 5 5.5733 1.61 1.052 74.9158-3 
4.3774 818 5 8 105 1.62 1.052 14.14353 
4.3874 7460 10 1 9.5111 1.61 ost wre 
4394 7430 4 3 2.515 1.61 1051 39.0145: 
4.3974 7.178 4 11 -L5'22 1.61 1.052 ало 
44074 6352 4 6 10723 1.60 1.053 18.3+ 1s 
4414 7519 5 4 3.017 1.60 1053 46.4+577 
4424 7436 11 2 100759 1.60 1053 156.217, 
44274 6788 2 2 1015 1.61 1053 т 
44374 71604 3 4  10'13 1.62 1.055 — 15.0223 
4444 7677 2 4 0051 1.63 1.055 Di 
4454 874 6 5 35119 1.64 1.055 45.51515 
44774 8167 7 5 45111 1.64 1.055 62.41372 
4.4974 7.997 2 9 —25t33 1.65 1055 —35.3t471 
45174 8.674 1 8 -30275 1.65 1055. 3897 
45314 9335 3 8 —1.0t3° 1.65 1055 = 
4.5474 8.765 1 5 —15725 1.65 1.054 —19.3738:8 
45514 8250 1 4 210128 1.65 1055 —13 71223 
4.5604 8390 2 1 157% 1.64 1055. ФӘ 
4.5774 8545 2 4 0075 1.64 1.055 (T 
4.5874 8162 4 3 2525 1.63 1055 34.74483 
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Definition of likelihood function 


In the maximum likelihood fit to the dressed cross sections of етет — ntr” Һе, the likelihood is defined as: 


17 62 
E(gP 5p) = [beep [T LÄS: ру), (1) 
4= 1 


j=1 


where 1/8 is the expected number of signal events, p; and p; are the parameters in the likelihood functions. L; and L; are the 
likelihood functions for the X Y Z and R-scan data samples, respectively. 


The likelihood functions for the X Y Z and R-scan data samples are defined differently due to the different statistics of the 
samples. For ће X Y Z data sample, where the statistics is large, the likelihood function is described by an asymmetric Gaussian 


; i ight 
function L; = С(иў®; mi, olf, o; 5"), 


(ш; 
і right 1 те і 
С(иЎї®; m; cleft о?! ) = - € 2(ciett)2 ue < т; 
à? PV wp left right TER ? 
у27(018 + оё") 
1 = WE- m)? 
2(oT ght 2 sig " 
a left, righty © ‘ B Z Mi; (2) 
V2n(of" +0; ) 


By scanning the number of Л. signal events in the fit to the yne invariant mass spectrum, the likelihood value as a function of 
expected signal events и; ° is obtained. The parameters of the Gaussian function are determined from a fit to the likelihood 
distribution. For the R-scan data sample, the likelihood is defined as: 


si si side S si si bk side bk bk 
Lut NS, Ме) =] PONG? uy + f - цо) РОМ, nj dp; 5, (3) 
0 


bkg 


where P(N; и) = oat uN e^" is the Poisson distribution, иў® and u; ° are the number of expected number of signal and 


background events, respectively. 


Correlation matrix from the fit to the cross section 


The correlation matrix of the fit parameters from the fit the to the dressed cross sections of e^ e^ — пел he 15 


Mi Г. (TeeB)1 ф M» T2 (TeeB)2 
М, 1.000 —0.420 0.104 0.136 0.184 —0.449 0.078 
Г! 1.000 0.191 0.077 0.478 0.096 0.079 
(TeeB)1 1.000 0.992 0.020 —0.102 0.985 
@ 1.000 —0.049 —0.119 0.993 
М» 1.000 —0.306  —0.113 
Гә 1.000 —0.017 
(Г.В): 1.000 


where М, Гу, апа (Teeb); (j = 1, 2) are the mass, the total width, and the product of the electronic partial and the branching 
fraction to ттлт he for the two resonances, respectively; ф is the relative phase between two Breit-Wigner functions. 


